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摘 要 : 为 探索 快速 提取 典型 绿洲 棉田 土壤 盐分 


新 疆 乌鲁木齐 ”830049; 3. 西北 农林 科技 大 学 水 利 与 建筑 工程 学 院 ,陕西 杨凌 712100) 


的 有 效 方法 ,获取 区 域 尺 度 的 土壤 盐 渍 化 特征 及 空 


间 分 布 ,进而 为 土壤 盐 渍 化 防治 提供 参考 。 以 新 疆 兵 团 农 二 师 31 团 为 研究 区 域 ,2019、2021 年 春季 
Landsat 8 OLI 多 光谱 影像 和 野外 实测 土壤 食盐 量 为 数据 源 , 将 波段 组 .光谱 指 数组 和 全 变量 组 作为 
模型 输入 变量 组 ,采用 多 元 逐步 回归 (Multiple stepwise regression, MSR)、 偏 最 小 二 乘 回归 (Partial 
least squares regression, PLSR ) 、 极 限 学 习 机 (Extreme learning machine, ELM)、 支 持 向 量 机 (Support 
vector machine, SVM) 和 BP 神经 网 络 (Back propagation neural network, BPNN ) 构建 基于 3 个 输入 变 


量 组 的 土壤 盐分 还 感 反 演 模 型 ,探究 输入 变量 和 建 模 方 法 对 模型 精度 的 影响 效果 ,通过 对 比 确定 


春季 土壤 盐分 最 优 反 演 模型 ,定量 反 演 地 表土 壤 含 盐 量 。 结 果 表 明 :(1) 研究 区 主要 为 非 盐 化 土 和 
轻 度 盐 化 土 , 总 样本 变异 系数 为 0.67, 呈 中 等 变异 性 ;光谱 反射 率 与 土壤 盐 渍 化 程度 的 关系 表现 为 
土壤 盐 渍 化 越 重 ,光谱 反射 率 越 高 。(2) 海岸 波段 (b1)、 蓝 波段 (b2)、 绿 波段 (b3)、 红 波段 (b4) 和 盐 
分 指数 (SI1 .SI2 .SI3 、SI4、S3、S4、S5) 均 通过 显著 性 检验 P<0.01, 相 关系 数 均 达到 0.4 以 上 。(3) 所 有 
模型 中 ,基于 全 变量 组 建立 的 BPNN 反 演 模型 精度 最 高 , 建 模 集 尼 为 0.705; 验 证 集 尼 为 0.3556。(4) 


由 反 演 结果 可 知 ,2019、2021 年 春季 耕作 区 土壤 


主要 为 非 盐 化 土 , 分 别 占 耕作 区 总 面积 的 $3.3$% 和 


64.62% ,其 次 为 轻 度 盐 化 土 , 分 别 占 44.31% 和 35.17% ;2021 年 土壤 盐 渍 化 程度 较 2019 年 有 所 减轻 。 


关 键 词 : 多 光谱 遥感 反 演 ; 土壤 盐分 ;光谱 反射 率 ; 
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目前 ,土壤 盐 渍 化 已 经 成 为 全 球 性 生态 环境 问 
题 ,严重 制约 着 干旱 区 灌溉 农业 的 可 持续 发 展 一 。 
新 疆 作 为 我 国 重要 的 农垦 区 , 盐 渍 土 面积 约 占 耕 地 
总 面积 的 1/3 , 盐 渍 化 问题 尤为 严重 ”。 因 此 , Ha 
快速 准确 获取 土壤 含 盐 量 变化 的 监测 方法 ,以 支持 
灌区 灌溉 制度 制定 、 土 壤 盐 涡 化 防治 .盐碱地 综合 
利用 土地 管理 等 工作 ”。 传 统 的 实地 采样 调查 方 
法 不 仅 耗费 大 量 的 人 力 ,物力 ,而 且 研 究 太 度 较 小 ， 
获取 的 多 为 点 状 信息 ““。 遥 感 技术 能 够 大 面积 
快速 获取 同一 地 区 不 同时 间 点 的 地 物 信 息 , 有 助 
于 开展 长 时 间 序 列 的 土壤 盐 溃 化 研究 工作 ,加 之 
宏观 .综合 、 及 时 动态、 低 成 本 等 特点 ,被 广泛 应 
用 于 土壤 盐 渍 化 的 实时 监测 和 土壤 盐分 的 定量 
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估算 等 1。 

随 着 遥感 技术 的 愈加 成 熟 , 学 者 们 开始 深入 探 
索 提高 建 模 精度 的 各 类 方法 ,其 中 ,光谱 参量 的 构 
建 和 建 模 方法 的 选取 有 助 于 提高 遥感 反 演 精度 ”。 
许多 学 者 采用 敏感 波段 组 合 ,特征 光谱 指数 构建 和 
光谱 数学 变换 等 方法 创建 模型 输入 变量 ”1。 周 晓 
红 等 ”以 匡 比 湖 湿 地 自然 保护 区 为 试验 区 ,基于 
Landsat 8 OLI 多 光谱 遥感 影像 ,构建 了 多 种 植被 指 
数 和 盐分 指数 ,在 众多 光谱 指数 中 ,增强 型 植被 指 
数 与 土壤 盐分 的 相关 性 最 高 ,传统 型 植被 指数 次 
之 ,盐分 指数 最 低 ; 张 晓 光 等 ”对 光谱 数据 分 别 进 
行 了 对 数 .倒数 .一 阶 微分 、 二 阶 微分 等 数学 变换 ， 
以 增强 光谱 与 盐分 之 间 的 相关 性 ,进而 达到 了 提升 
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模型 反 演 精度 的 效果 ;Zhang 等 “基于 MODIS 卫星 
影像 构建 了 光谱 指数 , 反 演 了 黄河 三 角 洲 地 区 的 土 
壤 含 盐 量 ; Alexakis 等 ”针对 希腊 克 里 特 岛 土壤 盐 
渍 化 问题 ,基于 WV2 All Landsat 8 OLI 计 算 盐 分 指 
数 ,结果 表明 ,光谱 指数 S5 对 土壤 盐分 最 为 敏感 ,可 
为 该 地 区 地 表土 壤 含 盐 量 的 监测 提供 新 方案 。 以 
上 研究 表明 ,筛选 敏感 光谱 变量 可 以 有 效 提 升 模型 
反 演 效果 ”。 除 了 自 变量 的 优化 , 建 模 方 法 也 经 历 
了 一 段 发 展 历程 ,从 传统 的 线性 回归 发 展 到 具有 非 
线性 拟 合 能 力 的 机 融 学 习 算法 。 人 研究 发 现 , 当 模型 
输入 变量 个 数 较 多 . 且 与 土壤 盐分 的 相关 性 较 低 
时 ,机 器 学 习 的 优势 得 到 了 更 加 充分 的 体现 ”。 通 
过 对 比 不 同 模型 的 建 模 和 验证 效果 ,确定 了 最 适合 
本 研究 区 土壤 盐分 的 建 模 方 法 和 反 演 模型 ,进而 提 
升 了 土壤 盐分 遥感 反 演 精度 。 冯 雪 力 等 “对 比 了 
传统 回归 方法 和 机 器 学 习 算 法 的 建 模 精 度 ,发 现 BP 
神经 网 络 (Back propagation neural network, BPNN) 
模型 预测 效果 最 优 ,其 稳定 性 和 预测 精度 均 优 于 经 
验 回 归 模 型 ; 刘 恩 等 % 基于 Landsat 8 0LI 多 光谱 影 
像 ,分 别 采用 多 元 线性 回归 和 BP 神经 网 络 算法 构建 
了 土壤 盐分 反 演 模型 ,发 现 影像 光谱 反射 率 和 土壤 
含 盐 量 之 间 不 是 单纯 的 线性 关系 ,因而 多 元 回归 模 
型 的 预测 效果 具有 很 大 的 波动 性 ,而 BPNN 具有 极 
强 的 非 线性 拟 合 能 力 ,因而 机 器 学 习 模 型 反 演 精度 
更 高 ; 张 智 揪 等 “以 内 蒙古 河套 灌区 为 研究 区 域 ， 
采用 多 种 方法 进行 建 模 , 发现 支持 向 量 机 (Support 
vector machine, SVM ) 模 型 精度 最 高 ,分 位 数 回归 次 
之 , 偏 最 小 二 乘 回归 (Partial least squares regression, 
PLSR ) 最 低 。 以 上 研究 都 说 明了 机 带 学 习 算 法 在 土 
SHEER) ER ST EP BEE 

土壤 盐 溃 化 相关 人 研究 虽然 取得 了 一 系列 较为 
理想 的 成 果 , 但 现 有 的 遥感 反 演 研究 多 是 基于 某 一 
特定 时 期 的 遥感 影像 展开 , 仅 能 反映 当地 该 时 期 的 
土壤 盐 渍 化 状况 ,模型 的 推广 性 不 强 。 而 且 , 盐 注 
化 程度 受 时 空 变化 影响 ,不 同位 置 的 土壤 含 盐 量 存 
在 差异 性 。 因 此 ,探讨 某 一 地 区 不 同年 份 的 土壤 盐 
澳 化 状况 极为 重要 一 ”。 鉴 于 此 ,本 文 以 新 疆 兵 财 
农 二 师 31 团 为 研究 区 域 , 基 于 2019、2021 年 的 
Landsat 8 OLI 卫星 遥感 影像 和 相应 时 期 地 表土 壤 含 
盐 量 ,利用 多 元 逐步 回归 (Multiple stepwise regres- 
sion, MSR) .PLSR 极限 学 习 机 (Extreme learning ma- 
chine, ELM) SVM 和 BPNN 构建 土壤 盐分 遥感 反 演 


模型 ,实现 棉田 土壤 盐分 的 定量 反 演 ,以 期 为 当地 
土壤 表层 含 盐 量 的 实时 监测 提供 参考 。 


1 材料 与 方法 


1.1 研究 区 概况 

新 疆 生 产 建 设 兵 团 农 二 师 31 团 位 于 天 山南 芳 
塔里木 盆地 东北 边缘 ,塔里木 河 与 孔 涛 河 两 河 下 游 
的 冲积 平原 上 ,地 处 巴 音 郭 楞 蒙 古 自治 州 尉犁 县 境 
内 ,西南 临 塔克拉玛干 沙漠 ,东北 与 库 姆 塔 格 沙漠 接 
HE , Hb SHEA fi A 85°24’ ~88°30’ E, 39°30! ~42°20'N. 
灌区 土壤 类 型 主要 为 盐 土 .碱土 和 风沙 土 。 团 场 深 
处 欧 亚 内 陆 腹地 , 属 北 温带 大 陆 性 荒漠 干旱 气候 ， 
降水 量 年 际 变化 大 ,多 年 平均 降水 量 53.3~62.7 mm, 
多 年 平均 莹 发 量 2273~2788 mm。 现 有 灌溉 面积 
7.5X10° hm? ,灌溉 方式 以 膜 下 滴灌 为 主 , 灌 溉 用 水 来 
自 塔里木 河 与 恰 拉 水 库 。 生 育 期 灌溉 量 2700~3450 
mhm”, 自 每 年 5 月 开始 ; 冬 灌 量 约 为 3600 m+ hm”, 
自 每 年 11 月 开始 ,使 用 大 水 漫灌 。 地 下 水 位 保持 在 
1~3 m 之 间 , 地 下 水 矿 化 度 2~7 gL'; 塔 里 木 河 与 水 
库 矿 化 度 0.47~1.53 gL'; 灌 溉 排水 矿 化 度 10.58~ 
18.36 g- L'o 
1.2 野外 数据 采集 与 处 理 

野外 采样 选择 在 棉花 播种 前 (2019 年 3 月 20 一 
28 H 2021 £4 H 71—12 日 ) 进 行 ,该 时 段 为 裸 土 期 ， 
受 植被 影响 较 小 。 实 地 勘察 后 ,按照 从 水 库 引水 处 
沿 干 渠 走 向 布设 采样 点 ,同时 考虑 交通 条 件 、 分 布 
位 置 .土壤 类 型 等 因素 。2019 年 和 2021 FRAKE 
的 样 点 个 数 分 别 为 33 个 和 33 个 (图 1)。 
1.2.1 表层 土壤 合 不 量 测定 考虑 到 土壤 水 盐 的 变 
异性 较 强 ,从 预 设 点 向 外 辐射 约 10 m 的 范围 内 再 选 
取 2 个 点 ,利用 3 点 取样 法 采集 地 表 0~20 cm 深度 的 
土壤 ,使 用 四 分 法 将 3 处 土壤 混合 ,作为 该 编号 位 置 
上 的 一 个 样品 ,随即 将 样品 装 入 已 编号 的 自封 袋 
中 ,手持 GPS 仪 测定 每 个 样 点 的 经 纬度 ,拍照 记录 
各 样 点 处 的 盐 涡 化 程度 .作物 类 型 .植被 覆盖 度 等 
信息 。 将 土 样 带 回 实验 室 后 , 挑 出 杂 物 ,风干 ,人 研 
磨 , 过 2 mm FL 425i, PRL 18 g 土 样 ,在 锥 形 瓶 中 按 土 
水 质量 比 为 1:5 配置 土壤 溶液 ,用 玻璃 棱 搅 拌 均匀 
后 固定 在 振荡 器 上 震荡 10 min ,使 盐分 充分 溶解 , 静 
置 后 使 用 电导 率 仪 测定 土壤 悬 浊 液 的 电导 率 值 
(EC, pS-em") ,根据 经 验 公 式 :y=0.0051x-0.5241( 决 
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图 例 (a) 研究 区 位 置 (b) 2019 年 采样 点 
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(c) 2021 年 采样 点 
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图 1 研究 区 位 置 及 土壤 样 点 分 布 示意 图 


Fig. 1 Sketch maps of the study area location and soil sample distribution 


E FRA R’=0.9534 ) E ERE ER (e+ ko)", 
1.2.2 表层 土壤 光谱 曲线 采集 与 处 理 ”选择 晴 衣 、 
无 云 .无 风 的 日 子 ,在 12:00 一 16:00 之 间 进 行 光 谱 测 
定 。 使 用 美国 PP Systems 公司 生产 的 UniSpec-SC 光 
谱 分 析 仪 (波段 范围 310~1130 nm ) 测 定 土 壤 光 谱 反 
射 率 。 样 点 观测 前 需要 进行 标准 白板 扫描 和 暗 扫 
措 操 作 ,光谱 采集 时 ,探头 设置 在 垂直 距离 地 面 50 
cm 处 ,从 而 缩小 仪器 视 场 角 ,以 减 小 地 面 背 景 及 周 
围 其 他 地 物 对 土壤 光谱 的 影响 。 对 不 同 盐 渍 化 等 
级 的 土壤 测定 多 条 光谱 曲线 ,每 个 等 级 测定 10 次 ， 
取 平 均值 后 作为 每 个 盐 渍 化 等 级 土壤 的 实测 高 光 
谱 数 据 。 土 壤 光 谱 采 集 完 毕 后 ,使 用 Multispec 5.1.5 
软件 对 数据 进行 处 理 , 得 到 310~1130 nm 波段 范围 
的 土壤 光谱 反射 率 ,在 Origin 2018 软件 中 绘制 土壤 
实测 高 光谱 反射 率 曲线 。 
1.3 遥感 影像 数据 的 获取 与 处 理 

多 光谱 数据 采用 Landsat 8 OLI 卫星 遥感 影像 ， 
数据 来 源 于 美国 地 质 调查 局 (United States Geologi- 
cal Survey, USGS) 官 网 (https://earthexplorer.usgs. 
gov!) ,行列 号 143/031, 空 间 分 辨 率 30 m, 重 访 周 期 
16 do 根据 与 采样 时 间 相 近 且 云 量 较 小 的 原则 ,本 
文选 取 2019 年 3 月 22 日 和 2021 年 3 月 20 日 的 影像 
作为 本 次 试验 的 多 光谱 数据 源 。 使 用 ENVI 5.3 软 
件 对 影像 数据 进行 辐射 定 标 .大气 校正 和 裁剪 等 预 
处 理 操作 ,实现 遥感 数据 由 数字 量化 值 到 辐射 亮度 
值 再 到 地 表 反 射 率 的 逐步 转换 ,再 通过 ArcGIS 10.2 
提取 各 采样 点 对 应 的 影像 多 光谱 反射 率 ,在 Excel 
2016 中 计算 各 类 型 光谱 指数 。 
14 光谱 参量 构建 

本 文选 取 了 7 个 波段 和 12 个 盐分 指数 作为 土 


壤 盐 分 反 演 模型 的 备 选 和 输入 变量 ,光谱 参量 及 计算 
公式 如 表 1 所 示 ,包括 海岸 波段 (b1) 、 蓝 波段 (b2)、 
绿 波段 (b3) 、 红 波段 (p4)、 近 红外 波段 (b5) .短波 红 
外 1 波段 (b6) .短波 红外 2 波段 (bp7) .盐分 指数 [ 归 
一 化 盐分 指数 (NDSI) 、SI-T、SIl、SI2、SI3、SI4、S1、 
S283 .S4,S5,S6]. 

R1 光谱 参量 及 相关 计算 公式 


Tab. 1 Spectral parameters and related formulas 


光谱 参量 计算 公式 参考 文献 
波段 bl CA 16] 
b2 B 16] 
b3 G 16] 
b4 R 16] 
b5 NIR 16] 
b6 SWIR1 16] 
b7 SWIR2 16] 
盐分 指数 NDSI RAR 14] 
SI-T 100 x (R- NIR) 14] 
su GXR 2] 
SI2 JG +R + NIR® 2] 
SI3 G+R 2] 
SI4 (B-G) +(G-R) [10] 
Sl BIR 2] 
S2 (B-R)/(B+R) 2] 
S3 (GX R)/B 2] 
S4 BxR 2] 
S5 (BxR)/G 2] 
S6 (Rx NIR)/G 2] 


IÈ: CA B.G, R.NIR SWIR1、SWIR2 分 别 为 海岸 波段 (bl)、 蓝 波段 
(b2).\ 绿 波段 (b3) 、 红 波段 (b4)、 近 红外 波段 (b5)、 短 波 红 外 1 波段 
(b6) ,短波 红外 2 波段 (b7) 的 光谱 反射 率 (%);NDSI 为 归 一 化 盐分 
指数 ;SI-T\SI1 SI2 SI3 .SI4 、S1 、S2 83 、S4 、S5 、S6 均 为 盐分 指数 。 
下 同 。 
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15 遥感 反 演 模型 构建 

考虑 到 土壤 盐分 受 水 分 的 影响 ,不 同时 期 的 灌 
溉 气温 和 降水 等 因素 会 影响 土壤 水 分 和 盐分 含 
量 ,2 期 土壤 含 盐 量 分 布 的 变异 性 存在 差异 。 别 除 
异常 值 后 ,将 2019 年 和 2021 年 春季 样本 数据 混合 ， 
按照 土壤 含 盐 量 由 小 到 大 的 顺序 排列 分 组 ,根据 建 
模 集 与 验证 集 2:1 的 比例 进行 等 间隔 取样 中 ,最 终 
确定 了 建 模 集 样本 51 个 ,验证 集 样本 26 个 。 对 建 
模 集 影像 的 波段 盐分 指数 与 土壤 含 盐 量 进行 皮尔 
HAH SELES AT (Pearson correlation coefficient, PCC) , 
HL ay PTET TAPE DET SS ae EE TT E RIESE 
谱 参 量 分 为 波段 组 .光谱 指数 组 和 全 变量 组 ,并 作 
为 输入 变量 组 参与 建 模 , 建 模 方 法 选用 MSR PLSR、 
ELM SVM 和 BPNN ,共计 15 个 土壤 盐分 反 演 模型 。 
建 模 方法 如 下 : 

(1) 多 元 逐步 回归 (MSR ) 

MSR 在 多 元 线性 回归 的 基础 上 进行 改进 , 它 能 
进一步 优化 公式 ,确定 各 个 自 变 量 对 于 因 变 量 预 测 
值 的 比重 ,剔除 不 重要 的 变量 “ ,是 土壤 盐分 遥感 
反 演 最 常用 的 回归 方法 之 一 下 。 本 人 研究 的 MSR 建 
模 与 分 析 基 于 IBM SPSS Statistics 25 完成 。 

(2) 偏 最 小 二 乘 回归 (PLSR ) 

PLSR 是 一 种 多 因 变 量 对 多 自 变量 的 回归 建 模 
方法 ,能 够 在 一 定 程 度 上 消除 参量 之 间 的 多 重 共 线 
性 ,通过 不 断 地 迭代 ,最 终 模型 将 包含 原 有 的 所 有 
自 变量 3。 本 文采 用 The Unscrambler X 10.4 软件 
进行 建 模 与 分 析 。 

(3) 极限 学 习 机 (ELM ) 

ELM 是 一 种 单 隐 层 前 馈 型 神经 网 络 , 具 有 泛 化 
能 力 强 、 人 工 干 预 小 、 学 习 速 度 快 等 优势 *, 其 特点 
是 输入 层 与 隐 含 层 之 间 的 连接 权 值 和 隐 含 层 神经 
元 的 浆 值 随机 生成 所 。 本 研究 基于 MATLAB 2018 
平台 进行 ELM 模型 构建 ,通过 反复 试验 ,最 终 确定 
了 波段 组 的 神经 网 络 模型 结构 (输入 层 神 经 元 个 
数 - 隐 含 层 神经 元 个 数 -输出 层 神经 元 个 数 ) 为 4- 
11-1, 光 谱 指 数组 的 结构 为 7-11-1, 全 变量 组 为 
11-7-1。 

(4) 支持 向 量 机 (SVM ) 

SVM 的 原理 是 建立 一 个 最 优 分 类 超 平面 ,使 其 
分 类 间隔 最 大 化 ”。 在 进行 非 线 性 问题 的 求解 时 
具有 泛 化 能 力 强 .全 局 最 优 、 结 构 化 风险 最 小 等 优 
势 ”…。 本 文采 用 LIBSVM 工 具 箱 在 MATLAB 2018 


软件 中 建立 SVM 模型 ,通过 测试 确定 了 最 适 于 本 实 
验 数 据 的 svmtrain 函数 参数 ,SVM 类 型 设置 为 e-SVR 
回归 , 核 函 数 选用 高 斯 径 向 基 核 函数 (RBF)。 惩 昼 
参数 (c) 和 核 函 数 参 数 (g) 根 据 样本 数据 进行 测试 ， 
最 终 确 定 了 波段 组 c 为 900,g 为 0.25; 光 谱 指 数组 c 
为 367,g 为 0.14; 全 变量 组 c 为 4,g 为 0.09。 

(5) BP 神经 网 络 (BPNN) 

BP 神经 网 络 (多 层 前 馈 式 误差 反 向 传播 神经 网 
络 ) 是 一 种 监督 学 习 模 型 , 它 通过 对 代表 性 样本 进 
行 学 习 训 练 ,从 而 掌握 研究 对 象 的 本 质 特 性 ,具有 
很 强 的 自 组 织 和 上 自 适应 能 力 ,能 够 模拟 任意 的 非 线 
性 输入 输出 关系 人”。 本 研究 基于 MATLAB 2018 平 
台 , 进 行 BPNN 的 结构 设计 。 将 建 模 集 的 波段 组 、 光 
谱 指 数组 和 全 变量 组 分 别 作为 模型 输入 因子 ,土壤 
盐分 作为 模型 输出 因子 ,建立 土壤 盐分 BPNN 反 演 
模型 。 经 过 反复 训练 ,最 终 确定 了 最 适 于 本 实验 数 
据 的 BPNN 模 型 参数 ,训练 函数 选用 poslin, 隐 含 层 
传递 隐 数 选用 purelin ,输出 层 函 数 选用 trainlm ,模型 
最 大 迭代 次 数 1000 次 ,训练 目标 最 小 误差 0.00001， 
波段 组 的 神经 网 络 结构 为 4-9-1, 光 谱 指 数组 的 结 
构 为 7-9-1 ,全 变量 组 为 11-9-1。 
1.6 模型 评价 指标 

为 量化 模型 的 拟 合 能 力 ,评价 指标 采用 决定 系 
数 (Coefficient of determination, R°) 和 均 方 根 误差 
(Root mean square error, RMSE ) ”对 模型 的 建 模 和 
验证 效果 进行 综合 评价 。 尼 越 接 近 1, 表 示 模 型 拟 
合 程度 越 好 ;RMSE 越 小 ,表示 实测 值 和 预测 值 误差 
越 小 。 因 此 , 尼 越 大 ,RMSE 越 小 ,说 明 模 型 精度 越 
高 , 反 演 效果 越 好 。 


2 结果 与 分 析 


2.1 土壤 盐分 特征 

根据 新 疆土 壤 盐 碱 化 程度 等 级 分 类 标准 ,将 
土壤 分 为 非 盐 化 土 (0~3 g: kg ') 、 轻 度 盐 化 土 (3~6 
gkg"!) 和 中 度 盐 化 土 (6~10 g.kg0)3 个 等 级 。 土 壤 
含 盐 量 的 分 析 结 果 如 表 2 所 示 ,总 样本 统计 结果 显 
示 , 非 盐 化 土 . 轻 度 盐 化 土 .中 度 盐 化 土 所 占 比 例 分 
别 为 64.94% .27.27% .7.79% ,变异 系数 为 0.67 ,表明 
棉田 土壤 含 盐 量 呈 中 等 变异 性 。 
2.2 实测 高 光谱 特征 

如 图 2 所 示 ,土壤 高 光谱 波长 范围 为 310~1130 
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表 2 土壤 盐分 的 描述 性 统计 分 析 


Tab. 2 Descriptive statistical analysis on Soil salinity 


样本 数 /个 土壤 合 款 量 /gkg” 
数据 集 — - - - — — 标准 差 ”变异 系数 
总 样本 非 盐 化 土  ” 轻 度 盐 化 土 ”中 度 盐 化 土 最 小 值 最 大 值 平均 值 
总 样本 77 50 21 6 0.37 11.82 3.02 2.03 0.67 
建 模 集 51 33 14 4 0.42 11.82 3.05 2.07 0.68 
验证 集 26 17 7 2 0.37 9.64 2.97 1.95 0.66 
表 3 光谱 参量 与 土壤 盐分 的 相关 
非 盐 化 土 X AF = 5 RRD 5+8 性 
0.7 ---- 轻 度 盐 化 土 Tab.3 Correlation between spectral parameters 
0.6 =- PRERE and soil salinity 
ar 光谱 参量 相关 系数 (7) 
x 0. 
=e 波段 oI 0.591" 
e o 2 0.550" 
0.3 3 0.539" 
0.2 04 0.525" 
61 5 0.297 
6 0.196 
0.0 1 1 1 1 1 1 1 1 1 J 7 0.329 
200 300 400 500 600 700 800 900 1000 1100 1200 i 
wK 盐分 指数 NDSI 0.304 
波长 /nm 
SI-T 0.197 
图 2 不 同 盐 渍 化 等 级 土壤 光谱 反射 率 曲线 SI 0.535" 
Fig. 2 Spectral reflectance curves of soil with different SI2 0.452” 
salinization grades SB 0.533" 
SI4 0.416" 
nm。 地 表土 壤 含 盐 量 不 同 导 致 土壤 反射 率 存 在 差 Sl -0.162 
异 , 但 不 同 盐 渍 化 等 级 土壤 的 光谱 曲线 相对 平行 ， a ee 
> r ; S3 0.497" 
形态 相似 ,变化 趋势 基本 一 致 , 随 着 波长 的 增加 , E a fe 
壤 高 光谱 反射 率先 减 小 后 增 大 。 波 长 在 310-373 j one 
nm 之 间 ,光谱 反射 率 随 波长 的 增加 而 迅速 减 小 ; 波 s6 0.269 


长 在 373~1130 nm 之 间 ,反射 率 随 波长 的 增加 而 增 
大 ,其 中 在 373~1050 nm 之 间 ,曲线 增长 较 缓 , 当 大 
于 1050 nm 时 ,反射 率 迅 速 增 大 。 整 体 来 看 ,土壤 盐 
渍 化 越 重 ,光谱 反射 率 越 高 。 
2.3 特征 光谱 参量 优选 

根据 皮尔 逊 相关 系数 界 值 表 一 ,样本 数 (n) 等 
于 50 时 , 当 相关 系数 (7) 的 绝对 值 大 于 0.354 时 , 表 
示 在 0.01 水 平 上 显著 ,光谱 参量 与 土壤 含 盐 量 的 
Pearson 相关 系数 如 表 3 所 示 ,筛选 出 敏感 波段 (P< 
0.01) 和 敏感 光谱 指数 (P<0.01)。b1、b2、b3、b4 波 有 段 
显著 性 检验 P<0.01, 且 7r 均 达到 0.5 以 上 ,与 土壤 含 
盐 量 的 相关 性 大 小 依次 为 bl1>b2>b3>b4。SIl、SI2、 
SI3 、SI4 S3 、S4、S5 盐分 指数 与 土壤 含 盐 量 均 呈 显著 
相关 (P<0.01) ,相关 性 大 小 依次 为 S4>S5>SI1>SI3> 
S3>SI2>SI4, 且 7 均 达 到 0.4 以 上 ,说 明 ST1 S12, S13, 


注 :*# 表 示 显 著 性 检验 P<0.01。 下 同 。 


SI4 S3 S4 .S5 在 一 定 程度 上 可 以 表征 土壤 含 盐 量 。 
2.4 土壤 盐分 反 演 模型 构建 

选择 与 土壤 盐分 显著 相关 (P<0.01) 的 bl、b2、 
b3、b4 波 段 作为 波段 组 ,选择 与 土壤 盐分 显著 相关 
(P<0.01) 的 SIl1、SI2、SI3、SI4、S3、S4、S5 盐分 指数 作 
为 光谱 指数 组 。 为 了 避免 模型 中 遗漏 重要 的 光谱 
参量 , 且 考 虑 不 同类 别 光 谱 参 量 对 模型 的 影响 ,将 
上 述 4 个 敏感 波段 和 7 个 敏感 光谱 指数 共同 作为 全 
变量 组 参与 建 模 。 
2.41 土壤 盐分 回归 模型 ”由 表 4 可 知 ,同一 变量 
组 作为 模型 输入 变量 下 ,MSR 模型 建 模 集 尼 均 高 于 
PLSR ,日 数值 均 大 于 0.340,RMSE 均 低 于 PLSR , 且 
数值 均 小 于 1.680, 说 明 MSR 建 模 效果 较 好 。 对 模 
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4 基于 不 同 变量 组 的 回归 模型 


Tab. 4 Regression models based on different groups of variables 


ARR 建 模 表达 式 建 模 集 (n=51) 验证 集 (n=26) 
方法 R RMSE R RMSE 
波段 组 MSR y=624.08b1-440.71b2-29.42 0.411 1.588 0.473 1.414 
PLSR  Y=40.75b1+35.51b2+25.32b3+19.94b4-21.19 0.314 1713 0.496 1.382 
光谱 MSR Y=-94.05S12+170.77S14+246.21S5-14.30 0.341 1.679 0.377 1.544 
指数 组 PLSR  Y=13.32S11+6.02S12+9.3S13+26.55S14+7.86S3+16S4+16.58S5-16.74 0.269 1.768 0.449 1.447 
全 变 MSR  Y=592.81b1-318.31b3-75.19S12+462.30S14-23.33 0.439 1.549 0.359 1.562 
量 组 PLSR  Y=15.29b1+13.33b2+9.5b3+7.48b4+8.53S11+3.86S12+5.96S13+ mea, “ana ee he 


17.01S14+5.0353+10.2554+10.6255- 18.73 


型 进行 精度 验证 时 ,发 现 波段 组 的 MSR 模型 和 
PLSR 模型 尼 相 近 , 均 大 于 0.472, RMSE 均 小 于 
1.415, 验证 效果 较 好 。 在 光谱 指数 组 和 全 变量 组 
中 ,PLSR 模型 验证 集 尼 均 高 于 MSR ,其 数值 均 大 于 
0.448 , RMSE 均 低 于 MSR ,其 数值 均 小 于 1.448 ,说 明 
PLSR 验证 效果 较 好 。 综 合 对 比 两 种 回归 方法 的 不 
同 变量 组 建 模 和 验证 效果 ,发 现 基 于 波段 组 建立 的 
MSR 模型 是 表 4 所 有 模型 中 最 优 的 土壤 盐分 反 演 
模型 。 

2.4.2 土壤 盐分 机 器 学 习 算 法 反 演 模型 ”如 表 5 所 
示 ,波段 组 和 光谱 指数 组 作为 模型 输入 变量 下 , 建 
模 集 的 各 类 机 器 学 习 模 型 让 均 大 于 0.616, 其 中 
SVM 模型 的 建 模 集 尼 均 达 到 0.768 以 上 ,RMSE 均 小 
于 0.998 , 建 模 精度 最 高 ,BPNN 模型 次 之 ,ELM 模型 
精度 最 低 。 对 模型 进行 精度 验证 时 ,发 现 波段 组 和 
光谱 指数 组 的 ELM SVM 模型 RMSE 较 大 。 在 全 变 
量 组 中 ,BPNN 模型 的 建 模 集 和 验证 集 尼 均 为 3 个 模 


注 :MSR 为 多 元 逐步 回归 ;PLSR 为 偏 最 小 二 乘 回 归 ; 尼 为 决定 系数 ;RMSE 为 均 方 根 误差 。 下 同 。 


别 为 1.133 和 1.409 ,模型 精度 最 高 ,ELM 模 型 次 之 ， 
SVM 模型 精度 最 低 。 因 此 ,在 3 种 机 器 学 习 算 法 中 ， 
BPNN 模型 最 适 于 本 实验 数据 。 综 合 对 比 3 种 算法 
的 不 同 变量 组 建 模 和 验证 效果 , 发现 基 于 全 变量 组 
建立 的 BPNN 模型 是 表 5 所 有 模型 中 最 优 的 土壤 盐 
分 反 演 模型 。 对 比 表 4 和 表 5 可 知 , 机 器 学 习 模 型 
精度 均 高 于 线性 回归 模型 ,因此 ,本 文 最 终 使 用 基 
于 全 变量 组 建立 的 BPNN 模型 来 反 演 31 团 春季 上 土 
壤 表 层 含 盐 量 。 
25 土壤 盐分 空间 分 布 反 演 

对 比 图 3 和 图 4 可 知 ,土壤 含 盐 量 为 6 oko VI 
上 的 土壤 均 分 布 于 灌区 中 部 的 盐碱地 带 上 , 含 盐 量 
为 0~3 g'kg 的 土壤 均 分 布 于 治 区 内 部 远离 蕊 地 的 
区 域 , 越 靠 近 荒 地 ,土壤 含 盐 量 越 高 。 由 此 可 知 , 土 
壤 盐 分 反 演 结果 与 Kriging 插 值 结果 基本 一 致 ,与 调 
查 情况 相符 ,表明 耕地 盐分 反 演 结果 的 准确 率 较 高 。 

如 表 6 所 示 ,2019 .2021 年 春季 耕作 区 土壤 主要 


型 中 的 最 高 值 ,分别 达到 了 0.705 和 0.556,RMSE 分 ”为 非 盐 化 土 ,分 别 占 耕 作 区 总 面积 的 55.55% 和 
表 5 基于 不 同 变量 组 的 机 器 学 习 模型 
Tab.S Machine learning models based on different variable groups 
ashe 建 模 方法 i 建 模 集 (n=51) 验证 集 (n=26) 
R RMSE R RMSE 
波段 组 ELM 0.617 1.279 0.533 1.998 
SVM 0.782 0.985 0.476 1.996 
BPNN 0.736 1.067 0.538 1.460 
光谱 指数 组 ELM 0.621 1.274 0.508 1.986 
SVM 0.769 0.997 0.409 2.001 
BPNN 0.730 1.080 0.459 1.567 
全 变量 组 ELM 0.540 1.403 0.490 1.565 
SVM 0.533 1.515 0.488 1.658 
BPNN 0.705 1.133 0.556 1.409 


注 :ELM、SVM 和 BPNN 分 别 为 极限 学 习 机 支持 向 量 机 和 BP 神经 网 络 。 
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Fig.3 Kriging interpolation distributions of soil salinity in tillage area 
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(b) 2021 年 春季 
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mm 3~4 «mm 6~10 mm m3~4 mm 6~10 mm 
图 4 耕作 区 土壤 含 盐 量 反 演 等 级 分 布 


Fig.4 Inversion grade distributions of soil salinity in tillage area 


表 6 不 同 盐分 层次 的 像 元 个 数 及 比例 


Tab.6 Number and proportion of pixels at different salinity levels 


土壤 类 型 ERER h/g kg" — “= 一 e 
像 元 个 数 面积 占 比 /% 像 元 个 数 面积 占 比 /% 
非 盐 化 土 (0~3 gkg'') 0~1 4508 0.35 14 0.00 
1~2 92995 TAT 366240 28.25 
2~3 622671 48.03 471548 36.37 
轻 度 盐 化 土 (3~6 oko) 3~4 428918 33.08 398125 30.71 
4~5 128548 9.92 47257 3.65 
5~6 16989 1.31 10556 0.81 
中 度 盐 化 土 (6~10 ekg") 6~10 1848 0.14 2737 0.21 


64.62% ,其 中 含 盐 量 为 2~3 g- ke AY EE mE 
最 大 ,1~2 eke 的 土壤 分 布 面积 次 之 ,0~1 gkg 的 
土壤 分 布 面积 最 小 。 轻 度 盐 化 土 面积 分 别 占 


44.31% 和 35.17% , 其 中 含 盐 量 为 3~4 g: kg 的 土壤 
分 布 面积 最 大 ,4~5 ko ' 的 土壤 分 布 面积 次 之 ,5~ 
6 eke 的 土壤 分 布 面积 最 小 。 对 比 2 a 土 壤 盐 渍 化 
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程度 ,发 现 2021 年 春季 非 盐 化 土 面积 相 比 2019 4E 
增加 了 9.07% , 轻 度 盐 化 土 面 积 缩小 了 9.14% ,中 度 
盐 化 土 面积 变化 较 小 。 

由 图 4 可 知 , 含 盐 量 小 于 2 gkg 的 土壤 主要 分 
布 于 恰 拉 水 库 附 近 、 耕 作 区 中 部 和 东南 部 。 赵 接近 
南北 两 侧 荒地 的 耕地 土壤 含 盐 量 越 高 ,东北 部 的 小 
块 耕地 盐分 偏 高 ,土壤 盐 渍 化 程度 主要 为 轻 度 盐 化 
土 。 中 度 盐 化 土 主要 分 布 于 灌区 中 部 的 盐 渍 化 地 
带 上 ,与 实际 影像 位 置 相 重 合 ,取得 了 较为 准确 的 
反 演 效果 。 此 外 ,东北 部 的 小 块 耕地 部 分 区 域 和 灌 
区 内 的 零散 地 块 上 也 有 少量 分 布 。2021 年 反 演 图 
中 绿色 更 深 , 即 含 盐 量 小 于 2 g.kge- 的 土壤 分 布 面积 
更 大 ,在 整个 灌区 都 有 体现 。2021 年 含 盐 量 大 于 4 
g .kg 的 土壤 分 布 面 积 相 比 于 2019 年 有 所 减 小 ,图 
上 表现 为 橙 红色 区 域 缩小 。 综 上 可 知 ,2021 年 土壤 
盐 渍 化 程度 较 2019 年 有 所 减轻 。 


3 讨论 


31 盐 渍 土 特征 分 析 及 光谱 参量 优选 

棉田 土壤 含 盐 量 呈 中 等 变异 性 ,说 明 其 易 受 气 
候 地形 .水文 `. 地 质 等 自然 因素 和 灌溉 、 翻 耕 .施肥 
等 人 为 因素 的 影响 3。 采样 时 间 为 2019 年 和 2021 
年 春季 ,2 a 自然 环境 差异 较 大 ,不 可 控 因 素 较 多 ， 
而 本 文 将 2a 春 季 数 据 混合 ,以 中 和 不 同年 份 的 样本 
数据 变异 性 ,再 进行 建 模 集 和 验证 集 划 分 ,更 具有 
合理 性 。 本 文 对 建 模 集 OLI 多 光谱 波段 .光谱 指数 
与 实测 土壤 盐分 进行 相关 性 分 析 , 发现 多 光谱 敏感 
波段 为 蓝 、 绿 、 红 波段 ,这 与 孙 亚 楠 等 人 研究 结果 一 
致 , 说 明 波段 为 430~760 nm 的 土壤 盐分 光谱 信息 量 
较 多 ,而 通过 波段 组 进行 特征 光谱 指数 的 构建 可 以 
挖 气 光 谱 与 盐分 之 间 的 隐 仿 信息 ,推动 土壤 盐分 遥 
感 定 量 反 演 研究 的 发 展 。 王 雪 梅 等 2 WHS! 
发 现 光谱 指数 与 土壤 盐分 之 间 存 在 一 定 的 相关 性 ， 
在 估算 应 用 方面 具有 很 大 贡献 。 本 研究 也 显示 , 
敏感 波段 与 敏感 光谱 指数 共同 作为 模型 输入 变量 
时 ,土壤 盐分 反 演 模型 的 建 模 和 验证 效果 最 优 。 
3.2 土壤 盐分 反 演 模型 建立 

本 文 基于 3 个 变量 组 ,使 用 5 种 方法 建立 了 15 
个 土壤 盐分 反 演 模型 ,发 现 选择 不 同 的 输入 变量 和 
建 模 方法 ,模型 效果 存在 差异 。 研 究 发 现 , 基 于 机 
器 学 习 模 型 的 精度 和 稳定 性 均 优 于 传统 线性 回归 ， 
这 是 由 于 土壤 盐 渍 化 发 生机 理 的 复杂 性 ,土壤 光谱 


特征 与 表层 土壤 含 盐 量 之 间 存 在 复杂 的 非 线性 关 
系 , 而 机 器 学 习 算 法 具有 极 强 的 非 线性 拟 合 能 力 和 
推广 能 力 ,适合 模拟 变量 内 部 机 制 的 复杂 关系 。 许 
多 学 者 研究 发 现 机 器 学 习 算 法 的 反 演 效果 优 于 传 
统 回归 模型 “2 ,与 本 文 的 研究 结果 一 致 。 冯 雪 力 
abs SE EAT PCR MSR 和 PLSR 建 模 时 ,发现 MSR 模 
型 模拟 全 盐 量 效果 最 优 ,本 研究 也 发 现 基 于 波段 组 
建立 的 MSR 模 型 是 6 个 线性 回归 模型 中 最 优 的 土壤 
盐分 反 演 模型 。 传 统 线性 回归 建 模 效 果 虽 次 于 机 
右 学 习 算 法 ,但 其 简便 易 懂 , 回 归 方 程 清楚 可 见 , 具 
有 一 定 的 运用 价值 ;机 器 学 习 模 型 可 以 极 大 的 提升 
模型 反 演 精度 ,也 具有 一 定 的 研究 价值 。 
3.3 土壤 盐分 时 空 变化 特征 

土壤 盐分 Kriging 捕 值 图 基于 各 采样 点 处 土壤 
盐分 的 实测 值 生成 ,数据 源 为 点 状 数据 ,数据 量 较 
小 ,成 果 图 较为 粗略 , 仅 能 大 致 反 应 土壤 盐分 的 分 
布 状况 。 而 遥感 技术 能 根据 土壤 盐分 反 演 模型 将 
影像 上 所 有 栅 格 点 处 的 光谱 反射 率 转化 为 土壤 含 
盐 量 ,由 此 达到 量化 效果 ,再 对 栅 格 点 数值 进行 范 
围 划 分 ,可 以 得 出 不 同 盐分 层次 的 棚 格 点 个 数 及 所 
占 比 例 。 对 比 2 a 土 壤 盐 渍 化 程度 ,发现 2021 EE 
坏 盐 渍 化 程度 较 2019 年 有 所 减弱 ,这 可 能 与 气温 、 
RRE .灌水 量 有 关 。 气 温 越 高 ,蒸发 作用 越 强 , 土 
壤 中 的 盐分 向 表层 运 移 ;灌水 量 越 大 ,对 土 体 淋 洗 
作用 越 强 ,土壤 表层 含 盐 量 越 低 。 而 且 当 地 非常 重 
视 盐 渍 土 的 治理 ,制定 了 合理 的 农业 灌溉 排水 制度 
和 措施 , 随 着 时 间 的 推移 ,取得 了 一 定 成 效 。 

已 有 研究 多 是 针对 某 一 年 份 的 盐 涡 土 进行 反 
演 ”*””, 但 由 于 每 年 环境 气候 多 变 , 且 栅 格 点 内 地 物 
类 型 复杂 ,使 得 不 同年 份 的 遥感 影像 数据 存在 较 大 
差异 。 本 研究 考虑 了 不 同年 份 自然 环境 的 差异 性 ， 
选用 2 a 春 季 数 据 参 与 建 模 和 验证 ,通过 不 断 调 整 模 
型 参数 ,建立 最 优 反 演 模型 ,对 春季 土壤 盐分 进行 
反 演 ,取得 了 和 良好 的 效果 ,但 本 次 研究 对 象 为 非 植 
被 覆盖 土 ,未 能 反 演 作物 生长 期 的 土壤 盐分 ,所 以 
下 一 步 可 以 尝试 研究 作物 光谱 与 土壤 含 盐 量 之 间 
的 关系 。 另 外 ,试验 区 自然 因素 复杂 多 变 , 限 制 了 
土壤 盐分 的 实时 监测 ,因此 反 演 方法 在 其 他 地 区 和 
其 他 季节 的 适用 性 还 有 待 进一步 研究 。 


4 结论 


(1) 棉田 土壤 主要 为 非 盐 化 土 和 轻 度 盐 化 土 ， 
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总 样本 变异 系数 为 0.67, 呈 中 等 变异 性 。 光 谱 反 射 
率 与 土壤 盐 渍 化 程度 的 关系 表现 为 土壤 盐 涡 化 越 
重 , 光 谱 反 射 率 越 高 。 

(2) 通过 土壤 盐分 与 OLI 多 光谱 影像 的 波段 、 
光谱 指数 进行 相关 性 分 析 , ÆI b1 b2, b3 .b4 波 段 
和 SIl 、SI2、SI3 、SI4、S3 、S4、S5 盐分 指数 均 通 过 显著 
性 检验 P<0.01, 相 关系 数 均 达到 0.4 以 上 。 

(3) 所 有 模型 中 ,基于 全 变量 组 建立 的 BPNN 
模型 反 演 精度 最 高 , 建 模 集 尼 为 0.705, RMSE 为 
1.133 ; 验证 集 尼 为 0.556,RMSE 为 1.409。 综 合 分 析 
建 模 和 验证 效果 ,选用 基于 全 变量 组 的 BPNN 模型 
来 反 演 春 季 棉 田 表层 土壤 含 盐 量 。 

(4) 由 反 演 结果 可 知 ,2019 .2021 年 春季 耕作 区 
ER EB HARE Et, ol ARE KO RY 
55.55% 和 64.62% , 其 次 为 轻 度 盐 化 土 , 分 别 占 
44.31% 和 35.17%。 越 接近 南北 两 侧 荒地 的 耕作 区 
土壤 含 盐 量 越 高 ,东北 部 的 小 块 耕地 盐分 偏 高 。 
2021 年 土壤 盐 渍 化 程度 较 2019 年 有 所 减轻 。 
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Inversion and validation of soil salinity based on multispectral 
remote sensing in typical oasis cotton field in spring 


LIU Xuhui'’, BAI Yungang’, CHAI Zhongping', ZHANG Jianghui’, 
DING Bangxin®, JIANG Zhu’ 


(1. College of Resources and Environment, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, China; 2. Xinjiang 
Research Institute of Water Resources and Hydropower, Urumqi 830049, Xinjiang, China; 3. College of Water Resources and 
Architectural Engineering, Northwest A & F University, Yangling 712100, Shaanxi, China) 


Abstract: This study aimed to explore an effective method for extracting soil salinity from cotton fields in oasis 
and determine the characteristics and spatial distribution of soil salinization on a regional scale to provide 
reference for soil salinization control. The 31* Regiment of the 2" Division of Xinjiang Production and 
Construction Corps was taken as the research area; Landsat 8 OLI multispectral images and field measurements 
of soil salinity in spring 2019 and 2021 were taken as data sources; and the band group, spectral index group, and 
total variable group were taken as the model input variable group. Multiple stepwise regression (MSR), partial 
least squares regression, extreme learning machine, support vector machine, and back propagation neural network 
(BPNN) were used to construct a remote sensing inversion model of soil salt based on the three input variable 
groups. Precision evaluation was conducted, and the effects of input variables and modeling methods on model 
accuracy were explored. The best inversion model of soil salt in spring was determined through comparison and 
quantitatively inverted surface soil salt content. Results showed that (1) the study area was mainly composed of 
non-salinized soil and slightly salinized soil, and the coefficient of variation of total samples was 0.67, implying 
moderate variability. The relationship between spectral reflectance and soil salinization is as follows: serious soil 
salinization leads to great spectral reflectance. (2) Significance tests were conducted on coastal band(b1), blue 
(b2), green(b3), and red(b4) and the salinity indices of SI1, SI2, SI3, SI4, S3, S4, and S5 (P<0.01). The obtained 
correlation coefficients were all above 0.4, which can represent soil salinity to a certain extent. (3) Among the six 
linear regression models, the MSR model established on the band group had the best inversion effect, and the 
BPNN inversion model established on the full variable group had the highest accuracy. (4) According to the 
inversion results, the soil in spring 2019 and 2021 was mainly nonsalinized soil accounting for 55.55% and 
64.62% of the total tillage area, respectively, followed by mild salinized soil accounting for 44.31% and 35.17%, 
respectively. Soil salinization in 2021 was reduced compared with that in 2019. 

Key words: multispectral remote sensing inversion; soil salinity; spectral reflectance; variable group; machine 


learning 


